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We present measurements of the branching fraction and time-dependent CP-violating asymmetries 
in B° -> K S K° S K° S decays based on 227 million Y(4S) -> BB decays collected with the BABAR 
detector at the PEP-II asymmetric-energy B Factory at SLAC. We obtain a branching fraction of 
(6.9±g;|±0.6) x lO" 6 , and CP asymmetries C = -0.34lg;|| ± 0.05 and S = — 0.71±g;|| ±0.04, where 
the first uncertainties are statistical and the second systematic. 

PACS numbers: 13.25.Hw, 13.25.-k, 14.40.Nd 



The amplitude of time-dependent CP violation (CPV) 
predicted for b — > ccs decays of neutral B mesons 
in the Standard Model (SM) is sin 2/3 where (3 = 
a,rg(—V c dV* b /V t dV t * b ) is the CP violating phase difference 
between mixing and decay amplitudes, with the ele- 
ments of the Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing matrix [1]. This prediction has been well tested 
at the B factories in recent years [2]. The SM also pre- 
dicts the amplitude of CPV in b — > sqq (q = d, s) decays, 
defined as sin2/3 e ff, to be approximately sin 2/3. However, 
since b — ► sqq decays are dominated by one-loop tran- 
sitions that can potentially accommodate large virtual 
particle masses, contributions from physics beyond the 
SM could invalidate this prediction, making these de- 
cays especially sensitive to new physics [3]. An active 
program has arisen to measure /3 e ff in as many b — ► sqq 
"penguin" modes as possible [4]. However, many of these 
final states are affected by additional SM physics contri- 
butions that obscure the measurement of /3 e ff [5], or arc 
not CP eigenstates. Two decays to CP cigenstates that 
have been noted as having small theoretical uncertainties 
in the measurement of /3 e g are B° — > 4>K® [6-8] (CP-odd) 
and B° -» K° s K° s K a s (CP-even) [9]. 

In this Letter we present a measurement of time- 
dependent CP- violating asymmetries in the decay B° — > 
K^KgKg, along with a measurement of the branching 
fraction (BF). Until recently the small branching frac- 
tion [10] and the absence of charged decay tracks orig- 
inating at the B° decay vertex have limited the ability 
to extract CP parameters from B° — ► K^K^K^. How- 
ever, techniques recently developed to deal with the re- 
construction of the B° decay vertex in B° — ► K®tt° have 
made this measurement possible [13]. 

The time-dependent CP asymmetry is obtained by 
measuring the proper-time difference At = tcp — t tag 
between a fully reconstructed decay B° — > K®K°K ° and 
the partially reconstructed tagging B meson (Stag). The 
asymmetry in the decay rate /+ (/_) when the tagging 
meson is a B° (B°) is given as 

e -|At|/r 

f ± (At) = (1) 

[1 ± 5" sin (Am d At) =F C cos (Am d At) ] , 

where the parameters C and S describe the amount of CP 
violation in decay and in the interference between decay 
with and without mixing, respectively. Neglecting CKM- 
suppresscd amplitudes, we expect S — — sin 2/3 and C = 
in the SM. 



The results presented here are based on 226.6 ± 2.5 
million T(4S) -> BB decays collected with the BABAR 
detector at the PEP-II asymmetric-energy e + e~ collider, 
located at the Stanford Linear Accelerator Center. The 
BABAR detector [11] provides charged-particle tracking 
through a combination of a five-layer double-sided sil- 
icon microstrip detector (SVT) and a 40-layer central 
drift chamber, both operating in a 1.5 T magnetic field 
Charged kaon and pion identification is achieved through 
measurements of particle energy-loss in the tracking sys- 
tem and Cherenkov cone angle in a detector of internally 
reflected Cherenkov light. A segmented CsI(Tl) electro- 
magnetic calorimeter provides photon detection and elec- 
tron identification. Finally, the instrumented flux return 
of the magnet allows discrimination of muons from pions. 

Candidates for B° — > K®K®K° are formed by com- 
bining three K® candidates in an event. We reconstruct 
K® — ► 7r + 7r~ candidates from pairs of oppositely charged 
tracks. The two-track combinations must form a ver- 
tex with a 7r + 7r~ invariant mass within 12MeV/c 2 (about 
4cr) of the nominal K° mass [12], a reconstructed flight 
distance between 0.2 and 40.0 cm from the beam spot 
in the plane transverse to the beam, and an angle be- 
tween the transverse flight direction and the transverse 
momentum vector of less than 200mrad. For each B 
candidate two nearly independent kinematic variables 
are computed, namely the beam-cnergy-substituted mass 
toes = \J{s/2 + Pi VB) 2 /Ef and the energy dif- 

ference AE = Eg - yfs/2. Here, (Ei,pi) is the four- 
vector of the initial e + e~ system, y/s is the center-of- 
mass energy, ps is the reconstructed momentum of the 
B° candidate, and E* B is its energy calculated in the 
e + e~ rest frame. For signal decays, the toes distribu- 
tion peaks near the B° mass with an rms deviation of 
about 2.5 MeV/c 2 and the AE distribution peaks near 
zero with an rms deviation of about 14 MeV. We select 
candidates within the window 5.22 < toes < 5-30 GeV/c 2 
and -120 < AE < 120 MeV, which includes the signal 
peak and a "sideband" region for background character- 
ization. 

The sample of B° — > K° S K° S K° S candidates is domi- 
nated by random K®K®K® combinations from e + e~ — > 
qq (q = u,d,s,c) fragmentation. Monte Carlo (MC) 
studies show that contributions from other B meson de- 
cays can be neglected. We exploit topological observables 
to discriminate the jet-like e + e~ — > qq events from the 
more uniformly distributed BB events. In the T(4S) rest 
frame we compute the angle 6? between the thrust axis 
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of the B° candidate and that of the remaining particles 
in the event. While | cos8 T \ is highly peaked near 1 for 
e + e~ — > qq events, it is nearly uniformly distributed for 
BB events. We require | cosO T \ < 0.9, eliminating ~ 68% 
of the background. In addition we use a Fisher discrim- 
inant variable (J 7 ), based on the momenta and angles of 
tracks in the event [13], in the maximum- likelihood fit 
described below. 

For the 1.4% of events with more than one candi- 
date we select the combination with the smallest x 2 — 
J2i( m i ~ m K°) 2 / a mi> wri ere TOj {m K o) is the measured 
(nominal K®) mass and a mi is the estimated uncertainty 
on the mass of the ith Kg candidate. We also remove 
all B° candidates that have a KgK® mass combination 
within 3d (45MeV/c 2 ) of the Xco or % c2 mass. While 
we expect few Xco and Xc2 ~* K®K® decays in our final 
sample, these are b — > ccs decays that would bias the 
CP- asymmetry measurement. 

We extract the results from unbinned maximum- 
likelihood fits to the kinematic, event shape (J-), and 
At variables. We maximize the logarithm of an extended 
likelihood function 

£ = e -(N S +NB) X JJ + NbV ^ t 

i 

where Vs and Vb are the probability density functions 
(PDFs) for signal (S) and continuum background (B), 
Nt is the total number of events, and N$ and Nb are the 
event yields to be determined from the fit. The product is 
over the selected events. The observables are sufficiently 
uncorrelated that we can construct the likelihoods as the 
products of one-dimensional PDFs. The PDFs for sig- 
nal are parameterized from signal MC events. For back- 
ground PDFs we determine the functional form from data 
in the sideband regions of the other observables where 
backgrounds dominate. We include these regions in the 
fitted sample and simultaneously extract the parameters 
of the background PDFs along with the fit results. 

For the branching fraction fit we use only the 
kinematic and event-shape variables (Vbf — 
V{m ES )'P(AE)'P(T)). There are two yields and six 
continuum PDF parameters floated in the fit. There 
are 721 K^K^K^ candidates that pass all the above 
criteria, and the fit to this data yields N$ = 88 ± 10 
events and Nb — 633 ± 26 events. Figure 1 shows the 
toes and AE distributions for these events with the 
results of the fit plotted as curves. As a check we also 
add a fit component for random combinatorial B back- 
ground, with PDF parameters determined from large 
MC samples. This fit finds 14 ± 11 candidates assigned 
to the B background. These candidates come from the 
continuum background; the signal yield changes by less 
than one candidate. A signal reconstruction efficiency 
of 5.6% is derived from a large MC sample in which 
the K® reconstruction efficiency is carefully matched 




m ES (GeV/c 2 ) AE (GeV) 



FIG. 1: Distribution of (a) 7TIbs and (b) AE for all events 
that pass the selections used for determining the branching 
fraction. The solid (dashed) curves are the PDF projections 
for the signal plus background (background only) from the 
fit. 



with that observed in large hadronic data samples. As- 
suming equal production rates of B°B and B + B~, we 
determine B(B° -> K°K°K°) = (6.9±g;| ± 0.6) x 10~ 6 . 

The largest systematic error (5%) for the branching 
fraction measurement comes from our uncertainty on the 
efficiency of reconstructing K® — > 7r + 7r~ decays. We de- 
termine uncertainties of 4% for the effect of the candi- 
date selection cuts and 5% for the parametrization of 
the PDFs used in the fit. The remaining uncertain- 
ties, including possible error in modeling the K^R^K® 
Dalitz plot distribution in determining the signal effi- 
ciency, combine to 2%. 

The CP-fit PDF for a given tagging category is V^p = 
'PBF'P c (At 7 crAt)£ c where e c is the tagging efficiency for 
tag category c. The total likelihood £ is the product of 
likelihoods for each tagging category, and the free param- 
eters are determined by maximizing the quantity In C. 
Along with the CPV asymmetries S and C, the fit ex- 
tracts e c for the background and other background pa- 
rameters. The background PDFs include parameters for 
the Af-resolution function 1Z and for asymmetries in the 
rate of B° versus B° tags. We extract 25 parameters 
from the CP fit. 

We use a neural network to determine the flavor of 
the -Btag meson from kinematic and particle-identification 
information [14]. Each event is assigned to one of six 
mutually exclusive tagging categories, designed to com- 
bine flavor tags with similar performance and At resolu- 
tion. We parameterize the performance of this algorithm 
with a data sample (Bn av ) of fully reconstructed B° — > 
Z)(*) _ 7r + / p + 1 a\ decays. The effective tagging efficiency 
obtained from this sample is Q = J2 C fC (l — 2w c ) 2 = 
0.305 ± 0.004, where e c and w c are the efficiencies and 
mistag probabilities, respectively, for events tagged in 
category c. 

We compute the proper-time difference At = {zcp — 
-Ztag)/7/?c using the known boost of the e + e~ system and 
the measured Az = zcp — zt ag , the difference of the re- 
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constructed decay vertex positions of the B° — > K®K®K® 
and Stag candidate along the boost direction (z). A de- 
scription of the inclusive reconstruction of the B tag vertex 
is given in Ref. [15]. For the B° -> K°K°K° decay, where 
no charged particles are present at the decay vertex, we 
constrain the B meson production vertex to the interac- 
tion point (IP) in the transverse plane using a geometric 
fit. The position and size of the interaction region are 
determined on a run-by-run basis from the spatial distri- 
bution of vertices from two-track events. The uncertainty 
on the IP position, which follows from the size of the in- 
teraction region, is about 150 /im horizontally and 4/zm 
vertically. The uncertainty on zcp, a convolution of the 
interaction region and the vertex of the B° — > KgKgK® 
decay, is about 75/im. The uncertainty on z tag is about 
200/^m and thus the uncertainty in Az is dominated by 
the uncertainty in the vertex of the tagging decay. The 
resulting resolution is comparable to that in B° — >■ J/tjj 
tf°[13]. 

Simulation studies show that the procedure we use to 
determine the vertex for a B° — ► K° s K a s K® decay pro- 
vides an unbiased estimate of zqp- The estimate of the 
At error in an event reflects the strong dependence of the 
zcp resolution on the number of SVT layers traversed by 
the K® decay daughters. However, essentially all events 
have at least one K® candidate for which both tracks have 
at least one hit in the inner three SVT layers (at radii 
from 3.2 cm to 5.4 cm). In this case the mean At resolu- 
tion is comparable to that in decays in which the vertex is 
directly reconstructed from charged particles originating 
at the B decay point [15]. For a small fraction (0.1%) of 
the signal events, at least one K® has tracks with hits in 
the outer two SVT layers (at radii 9.1 cm to 14.4 cm) but 
none of the three s have hits in the inner three layers. 
In this case the resolution is nearly two times worse but 
the event can still be used in the CP fit. Events with 
a At > 2.5 ps or |Ai| > 20 ps are excluded from the CP 
fit. 

The resolution function 1Z is parameterized as the sum 
of a 'core' and a 'tail' Gaussian distribution, each with a 
width and mean proportional to a At, and a third Gaus- 
sian with a mean of zero and a width fixed at 8ps [15]. 
We have verified with MC simulation that the parame- 
ters of 1Z for B° — > K®K2.K2, decays are similar to those 
obtained from the -Bflav sample. Therefore, we extract 
these parameters from a fit to the -Bfl av sample. We find 
that the At distribution of background candidates is well 
described by a delta function convolved with a resolution 
function having the same functional form as that for the 
signal. The parameters of the background function are 
determined in the fit. 

The fit including At and tagging information yields 
S = — 0-71±g:|| ± 0.04 and C = -0.34±g;|? ± 0.05. Fix- 
ing C = we obtain sin 2/3 = -S = 0.79±gJ§ ± 0.04. 
Figure 2 shows distributions of At for £?°-tagged 
and -B°-tagged events, and the asymmetry A(At) = 




-6 -4 -2 2 4 6 

A t (ps) 

FIG. 2: Distributions of At for events weighted with the sPlot 
technique for B ta g tagged as (a) B° or (b) B°, and (c) the 
asymmetry _4(At). The points are weighted data and the 
curves are the PDF projections. 

(N B o - Ng ) I (N B o + N So ) , obtained with the sPlot 
event weighting technique [16]. 

TABLE I: Systematic uncertainties on S and C. 

Resolution function 0.017 0.017 

Vertex reconstruction 0.020 0.022 

SVT alignment 0.015 0.008 

Background asymmetry 0.007 0.022 

Fit correlation 0.016 0.004 

Tag-side interference 0.008 0.015 

PDFs 0.025 0.026 

Tbtal 0.044 0.047 



Systematic uncertainties on the CP parameters are 
given in Table I. The systematic errors are evaluated with 
large samples of simulated Bflav and B° — > K®K®K® 
decays. We employ the difference in resolution func- 
tion parameters extracted from these samples to vary 
the resolution function parameters extracted from the 
Bflav sample in data. We also perform fits to the simu- 
lated B° — > K®K®K® signal with parameters obtained 
either from signal or i?fl av events to account for any 
potential bias due to the vertexing technique. Several 
SVT misalignment scenarios are applied to the simulated 
B° — > KgKgKg events to estimate detector effects. We 
consider large variations of the IP position and resolu- 
tion and find they have negligible impact. Asymmetries 
in the rate of B° versus B° tags in the background events, 
which are free parameters in the fit, are fixed to zero as 
a systematic uncertainty. The systematic error due to 
correlations in the fit variables is extracted from a fit to 
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a sample of randomly selected signal MC events added 
to background events from a parametrized MC. We al- 
low for the possible interference between the suppressed 
b — > ucd and the favored b — > cud amplitude for some 
tag-side B decays [17]. Finally, we include a systematic 
uncertainty to account for imperfect knowledge of the 
PDFs used in the fit. Most of the uncertainties on the 
PDFs arc statistical and some are associated with data 
and MC differences. 

In summary, we have measured the B° -» K°K°K° 
branching fraction and the time-dependent CPV asym- 
metries. The BF measurement is in good agreement with 
previous measurements [10]. The measurements of S and 
C are in good agreement with the SM expectation. 
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